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Abstract 
in 

■^■j- ■ We inspect the first-order electron-electron capture scenario for transfer ionization that has been 

o: 

recently formulated by Voitkiv et al. (Phys. Rev. A 86, 012709 (2012) and references therein). 
f ^ | Using the multichannel scattering theory for many-body systems with Coulomb interactions, we 

m ■ 

show that this scenario is just a part of the well-studied Oppenheimer-Brinkmann-Kramers ap- 
proximation. Accurate numerical calculations for the proton-helium transfer ionization reaction in 

X ■ 
j_j \ the first Born approximation exhibit no appreciable manifestation of the claimed mechanism. 
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I. INTRODUCTION 

Various capture processes, including single and double charge transfer, transfer excitation 
(TE), and transfer ionization (TI), have been attracting interest from researchers for nearly 
a hundred years. Recent H + + He experiments allow to distinguish not only the ground 
state of the residual ion He + , but also its excited and ionization states, as well as the 
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-|4|]. To our knowledge, the first 



ground and excited states of the final atomic hydrogen 

(classical) description of the capture mechanism was given by Thomas in 1927 J5|. Briefly, 
an electron, after being kicked by a fast proton, elastically scatters from the target nucleus. 
The scattered electron can move parallel to the proton with the same velocity. This is a 
classical description of the capture process. To transfer the electron into a hydrogen bound 
state, a gamma-quant must be emitted to carry away the released energy (the so-called 
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radiative capture |6|). This mechanism is referred to as Ne-Th.ora.as. 

Almost at the same time, a quantum mechanical explanation how a target electron can be 
captured into a bound state of a fast moving projectile (proton), was given by Oppenheimer, 
Brinkmann and Kramers (OBK) [7J]. In the OBK scenario the electron transfer proceeds via 
an overlap of initial and final wave functions of the projectile-target system. This so-called 
kinematical capture relies strongly on the radial and angular electron correlations in the 
target, if we consider TE and TI processes. It must be noted that in quantum mechanics 
the transition of an electron to the hydrogen bound state can be nonradiative, and the 
energy excess can be carried away by a third body that participates in the reaction. In 
1974, Shakeshaft J8j demonstrated connection between the correct classical and quantum 
mechanical second-order treatments of the iVe-Thomas process at high projectile velocities. 
And later, Briggs and Taulbjerg 9[ carried out a quantum treatment of the so-called ee- 
Thomas mechanism, where an electron after the first kick is scattered by another target 
electron, not by the nucleus. 

In quantum mechanics fast processes are usually treated in terms of Born approxima- 
tions. This framework is directly applicable in the case of two-body scattering, but requires 
additional careful considerations in the many-body case. The situation of particular im- 
portance is when the entrance channel of the reaction is different from its exit channel, for 
example, as it is in capture processes. Within the multichannel scattering theory, the OBK 
mechanism can be attributed to the first Born approximation (FBA), whereas the Ne- and 



ee-Thomas mechanisms can be described using the second Born approximation (SBA). Any 
Born approximation is a sum of matrix elements. Each of them corresponds to a particular 
interaction that enters a total perturbation potential. For example, the OBK matrix element 
is one of the three FBA terms. The SBA contains twelve different terms, and only two of 
them correspond to the Ne- and ee-Thomas mechanisms. 



Recently Voitkiv et al. published a series of papers |l0l-ll2j putting forth a new first-order 
capture mechanism that can be called electron-electron Auger (or ee- Auger) 12[. According 
to this mechanism, the electron undergoes a nonradiative transition from the atomic state 
to the bound state of the projectile, transferring the energy excess to the another atomic 
electron which is emitted from the atom. This scenario resembles a kind of Auger decay and 



to be contrasted with the first-order radiative capture [13| which is accompanied by emission 
of a photon instead of an electron. A clear signature of the ee- Auger mechanism, according 
to Voitkiv et al, is emission of the electron in the direction opposite to the projectile motion 
(in the rest frame of the atom). In this work, we examine the mechanism suggested by 
Voitkiv et al. on the basis of consistent multichannel scattering theory. We show that, in 



contrast to the claim of Ref. [10|, it is not new and previously undiscussed. Namely, it is 
just a part of the usual kinematic capture in the OBK approximation, and the correlated 
electron-electron emission is nothing else but a result of the electron-electron correlations in 
the target atom. 

The paper is organized as follows. Sec. [Ill is devoted to quantum scattering theory of the 
considered transfer ionization process. In Sec. IHIl effects of distortion beyond the first-order 
approximation are discussed. Numerical results are presented in Sec. [TV] Finally, Sec. IVl 
summarizes this work. Atomic units (a.u., H — e — m e = 1) are used throughout unless 
otherwise specified. 

II. MULTICHANNEL SCATTERING THEORY 

In this section we remind basic formulas of quantum scattering theory for many-body 
systems. More mathematical details, particularly for the case involving charged fragments, 
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one can find in the review article [14J. A set of relative momenta defining motion of n 
fragments colliding in the asymptotic channel is denoted by p a . In turn, the ket vector 
\cf) a ) stands for a product of bound (spectral) state wave functions, which define the channel 



a. Hence, the ket vector \(j) a ,p a ) is the eigenfunction of the asymptotic hamiltonian H a : 
(E a — H Q )\(j) a ,p a ) = 0. The total hamiltonian is H = H a + V a , where V a is a set of two-body 
interaction potentials, which we consider as perturbation, and they define the terms of the 
Born series: FBA, SBA, and so on. 

The amplitude of the transition from the channel a to the channel /3 can be presented 
using two forms. These are the post-form 

T Pa {E) = (^,p>|^ + ^(E-i/ + zO)- 1 K|0 a ,p Q ), (1) 

and the prior-form 

f Pa {E) = (<f>p,pp\V a + Vp{E-H + iO)- l V a \<f> a ,p a ). (2) 

It is straightforward to demonstrate that on the energy shell, where E = E a = Ep, one has 



(see, for instance, Ref. |l5J |) 

Tp a (E) = fp a (E). (3) 

Moreover, since V a = H — H a and Vp — H — Hp, the relation (3) holds true in the FBA 
case as well, that is, on the energy shell the FBA post- and prior-amplitudes coincide, 

{fo>Pp\Vfi\4> a ,Pa) = (<(>f3,Pp\V a \4>a,Pa)- 

In the case of charged fragments the scattering theory should be modified. This modi- 
fication is due to the fact that charged fragments continue to interact even at asymptoti- 
cally large distances from each other, and their trajectories are not straight lines anymore. 
This asymptotic interaction in the channel a is characterized by the so-called Sommerfeld 
(Coulomb) parameter 
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Here Z, is the total charge of the ith asymptotic fragment (ion, atom, and electron), and % 
the relative velocity of fragments % and j. At this point we introduce distorting potentials 
and distorted waves. There are many ways of doing this, but a channel distorted function 



Xa (Pa) must obey the asymptotic relation at t — > =Foo [16] 



\X*(Pa)) ~+ e-^^^^l^pa). (4) 



Here A a (p a ) is the Dollard phase, whose explicit form is not important for the following. A 
distorting potential (operator) must satisfy the only condition 

and the channel wave function ^f^{p a ) is given by the modified Lippmann-Schwinger repre- 
sentation (jj, Ixrjl 

l*a(Pa)> = IxtiPa)) + (E a ~ H ± iO)" 1 ^^ + [H, a , V±] ) 1 Q , p a ) . (5) 



All well known approximations (the Peterkop effective charges 18J, CDW 19|, |20|, and etc.) 
are related to Eq. (5). The only difference between them is the choice of the operator V^ 
Thus, we can present the post- and prior-amplitudes as (J3 7^ a, E — E a = E^ 



3 



T^E) = QpMWW + [#* V^D+ltf+fo)), ( 6 ) 

f Pa {E) = {V-p&WaYt + [H a , V+]|0 Q ,p Q ). (7) 

These relations reduce to Eqs. (1) and (2), respectively, when V^ = 1 or, equivalently, 
Va,/3 = 0. As previously, the post- and prior-amplitudes coincide. Another important re- 
mark is that, if we account exactly for a particular distorting potential in the channel a, 
for example, for the Coulomb potential between two fragments, we must subtract the cor- 
responding Sommerfeld parameter from the total channel parameter r] a (see, for instance, 
Ref. (Jl). 

Let us apply the above general formulas to the fast TI reaction H + + He — > H + He 2+ + e 



discussed in the papers of Voitkiv et al. [KH12J. The authors utilize the pos£-amplitude, 



which in the nonsymmetrized FBA can be written as 

T£ BA (E) = (4> pl , tpN 2 (k),p H \V m + V p2 + V 12 + V Np \%,p ). (8a) 

In Eq. (8a), electrons are labelled by "1" and "2", "p" labels the fast proton projectile, and 
"iV" the target nucleus. The wave function \<j) p i) is the bound (ground) state of atomic 
hydrogen, \(fN2{k)) the continuum state of the He + ion, |$o) the helium wave function, po 
the proton momentum, pn the hydrogen momentum, and k the momentum of the emitted 
electron. Further, Voitkiv et al. put forth their ee-Auger mechanism by leaving only the 
term V\ 2 in the matrix element (8a). At the same time, they provide no evidence that this 



specific term indeed dominates in the FBA amplitude (8a). Moreover, taking into account 
that this amplitude is equal to that in the prior-form 

f™ A (E) = (</> p i,<pJr 2 {k),PH\V pl + V p2 + V Np \$ ,p ), (8b) 



there is no ground for such a claim. It was shown many times (see, for instance, Ref. 23]) that 
in the FBA prior-amplitude (8b) even the OBK term, corresponding to V p \ and containing 
the ee- Auger mechanism (see below), is not dominant in some kinematical situations. In 
other words, all four terms in (8a) must be taken into account. From the equality of the 
FBA amplitudes (8a) and (8b) we find that 

{<PpUVn2^)^Ph\V p i\^q,Pq) = ((j> p l,(PN 2 ( k )iPH\ V Nl + Vi 2 \®0,Po}- (9) 

It follows that the OBK matrix element, occurring in the left-hand side, is equal to the 
sum of the two matrix elements generated by the potentials V/vi and V12, respectively. 
Thus, it appears that the ee-Auger mechanism claimed by Voitkiv et al. as new and so far 
unconsidered is, in fact, already included in the very well studied OBK scenario. 

Eqs. (8a) and (8b), are valid when the total Sommerfeld parameters in the initial and 
final channels, t] a and 77,3, are equal to zero. In the present case, it is clearly fulfilled in 
the initial channel (the He atom is neutral and does not interact with the charged proton 
at asymptotically large separation distances). But it is also fulfilled in the final channel, 
because we use in (8) the spectral Coulomb functions <^^i 2 {k) instead of plain waves and, as 
remarked above, the Sommerfeld parameter of the N2 subsystem does not contribute to the 
total one. 

Now we consider the Schrodinger equation for the discussed problem in real space. Below 
we shortly recapitulate results from Ref. [2J. Let us introduce spatial positions of electrons 
in the laboratory frame r 1; r 2 , and R is the position of the projectile proton. If we consider 
very small scattering angles (fractions of milliradians) and not too high projectile energies 
(within few MeV), then the target nucleus stays practically at rest in the laboratory frame 
during the collision and we can set its position as the origin of the coordinate system, i.e. 
rV ~ 0. We also introduce the position of electron 1 relative to the proton x = R — T\. 
After a simple paraxial transformation, ^/j(p H , k\ R, x, r 2 ) — » exp(ip H ■ _R)$7(/c; R, x, r 2 ), we 



obtain for the wave function |<&. ) the following equation: 
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\&2 + \-V N2 ) \<f>](k;R,x,f 2 ) 



(V p2 + V p n + V m + V12) ^7(fc; R, x, r 2 ) 



(10) 



with the final boundary condition |$7) — ^ |0 P i, V^Tv^^))- Since the proton mass is m = 1836, 
the term A R /2(m + 1) can be neglected in (10). Setting R z = v p t and R± = b, i.e., 
R(t) = Vpt + b, we obtain the effective "time- dependent" Schrodinger equation 
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12 2 1 
\R(t)-r 2 \ R(t) \R{t)-x\ \R(t)-x-r 2 \ 



In this connection, let us remind a well known relation 
independent scattering theory and the approach of Refs. 10 



$j(k;b,x,r 2 ;t) = 

<$>j(k;b,x,r 2 ;t). (11) 

between the present time- 
12j, where time-dependent 



formulation with classical treatment of the fast moving projectile is employed. In that 
approach, the heavy incoming ion (proton) is considered as a source of a classical field en- 
tering the time- dependent Schrodinger equation, and its trajectory in the laboratory frame 
is R(t) = b + Vpt, (b ■ v p ) = 0, where b is the impact parameter. The time-dependent 
Schrodinger equation is 

.0 



'at 



$*(*)) = [H + V pl (t) + V p2 (t) + V Np (t) + V m + V N2 + V 12 ]\$ ± (t)), 



(12) 



where Ho is the kinetic-energy part of the total hamiltonian H. The interaction potential 
of the proton with the helium target is V p He(t) = V p \(t) + V p2 (t) + VJv p (£) ~ 0(l/t 2 ), if 
t — > ±00. The initial and final conditions at t — > =Foo, respectively, are 



(fi,r 2 |<l> + (£ -)■ -00)) -¥ e 



-ie" e t 



[ri,r 2 \® ), 



(x,r 2 \$-(t -> +00)) -> e ^/2-^-^ 2 /2)* e -^-(f,f 2 |0 pl5V ,- 2 (fc)), 

where e^ e ~ —2.9037 is the helium ground-state energy, e H = —0.5 the ground energy of 
hydrogen atom, and x = R(t) — f\. The so-called translational factor e l ^ v p' 2 ' t ~ Vp '^ J9( arises 
as a result of Galilean transformation of the hydrogen state from the projectile-fixed to the 
laboratory (target-fixed) frame. 



Bearing in mind the initial condition, the asymptotic time-independent hamiltonian in 
the initial channel is taken as H°- s = H + Vni + Vj\T2 + V12, with Wi(t) = V P He(t) being 
perturbation. And in the final channel we have Hf s (t) = H + V p i(t) + Vjv 2 and ^/(t) = 
V P 2{t) + VJvp(t) + Vjvi + V\2- The transition prior- amplitude in the time-dependent theory is 

/oo 
dt($-(t)\W t (t)\®o)e- t£ ° H , (13a) 

-oo 

while the pos£-amplitude is given by 

/oo 
dt(cW*'«^ 1 ,^(ft|W>(t)|* + (t)>e*^ +fca ^, (136) 

-00 

where e l ^ Vp '^ is the translational factor (see above). Both forms, (13a) and (13b), are 

equivalent, provided $ ± (t) is the exact solution of Eq. (10) 15|, |20[. Note that the authors 



of Refs. 
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- 12| unjustifiably set Wj(t) = V12 in the amplitude (13b). 
Comparing to Eq. (11), one readily sees the equivalence of the time-dependent and 
stationary formulations of scattering theory in the paraxial approximation. Finally, the 
relation between the T/j and a/j amplitudes is given by 

T fi = v p f d 2 b e~^ rb a fl (b), 

where q= p H ~Po, q± = q-q z v p /v p , v p = {0, 0, v p }, and k = {k± cos <p k , k± sin <p k , k z }. The q z 
value is determined from energy conservation as q z = v p /2+Q/v p , where Q = e^ e —e^ — k 2 j2. 

III. DISTORTED WAVE APPROXIMATIONS 

From Eqs.(l), (2), (6), and (7) one can derive the higher Born terms as well as different 
versions of the distorted wave Born approximation (DWBA). They include, in particular, 
matrix elements which correspond to both iVe- and ee-Thomas mechanisms, where, however, 
all two-body potentials must be taken into account (see, for instance, different calculations 

n n 

presented in Ref. |2j). For example, in Ref. [2[ the eikonal approximation was obtained, 
which introduces in the prior-form FBA matrix element (8b) a distorting phase factor, 

\</> P i,(p^(k),p H ) -» e {i/v ^\(j) pl ,ip], 2 (k),p H ). 

The details concerning its derivation on the basis of Eq. (11) in the closure approximation 
can be found in Ref. 2j. The result for the phase factor reads 

£ = ln ( WM) -x\+vy (R(t) ~ x)T [v P \R(t) ~ r 2 \ + v P ■ (R(t) - r 2 ) ]' 
/_ \ [v p R{t) + v p ■ R{t)} 2 [v p \R(t) - x - f 2 \ + v p • (R{t) - x - r 2 )} 



We note that it is the asymptotic form of the product 



Ji/vp)6f _vA-_A-A- A- A- 

e —fivj— i\ p2 i\ pN i\ Nl i\ 12 , 



(14) 



with 
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Z\Z% 



cxp 



2v„ 



r 1 



.Z\Z 2 , 

Vn 



.Z\Z 2 
% , 1; 

Vr, 



-i(v p r re i + v p f re i) 



Each factor ^■z 1 z 2 m (l^) describes the distortion of interactions between different con- 
stituents of the two final compound systems, pi and N2. In some sense, it is a 4C model 
(in analogy with well known 3C and 6C models in the scattering theory [22|). 

The same procedure we can utilize in the case of the FBA matrix element in the post-iorm, 
replacing in (8a) 



Here again 



p {i/v p )8 i A+ _ A + A+ A + 



(15) 



and 



Z1Z2 



cxp 



7TZiZ 2 

2v r , 



r 1 



Z\Z 2 
% ) 1F1 



- Z\Z 2 _ _ _ 

-% , 1, 1i\VpT re i VpTrel) 



P / \ P 

Each factor A^ iZa in (15) describes the distortion due to interactions between the projectile 
proton and different constituents of the helium atom. This approximation is analogous to 
the 3C model. 

Note that in the matrix elements we return to the time-independent formulation by 
replacing v p t — > R z . Both distortion factors, A7 and A^", partially incorporate an effect 
of higher Born terms. However, they principally differ from those employed in the papers 
of Voitkiv et al.. The latter distortion factors clearly do not obey the correct boundary 
conditions A7 — )■ 1 and A^ — > 1 at t —$■ ±00, respectively, and lack justification from the 
multichannel scattering theory. 



IV. NUMERICAL ILLUSTRATION 



In Sec. |TT]it was shown that in Refs. 1CH12| the authors use the FBA matrix element in its 
post- form, leaving only the V12 contribution. Their account for distortion effects beyond the 
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FBA through the factors A,;j is, at least, questionable (see the above remarks). Moreover, 
in Ref. 12] the final distortion factor was neglected in calculations, A/ = 1, as being not so 
much significant. As remarked in Ref. 12j, without the initial distortion factor, i.e., when 
Aj = 1, the contribution of the ee- Auger mechanism calculated there becomes much larger, 
while that of ee-Thomas vanishes. In view of these remarks, we might expect that neglecting 
the distortion effects in calculations does not reduce the role of the ee- Auger mechanism. 
Thus, in the FBA calculations presented in this section we set Aj = A/ = 1. Further, we 

1, not 



use FBA in the prior- form (8b), pointing out once again that it is equivalent to the fu 
truncated version of the post-iorm (8a). The quantity that was studied in Refs. 10Ml2J is 
the double differential cross section (DDCS) 

dip k / d 2 q±\T fi \ 2 , (16) 



d 2 a 



2k 



dk i dk, 



(2tt)V 



which describes a 2D distribution of the momentum components of the emitted electron 
(k x = k±cos<p k , k y = k ± simp k ). 

Fig. 1 shows numerical results for the DDCS using the FBA amplitude in the prior- 
form (8b) with an accurate, highly correlated helium function of Mitroy et al. 24] . The 
kinematical situation is the same as that of Fig. 1 in Ref. 12]. However, while there is 
a common feature such as maximum located at k z = 0, we find no maximum located at 
negative k z values (approximately, at k z « —3.0), in contrast to the results presented in 
Fig. 1 of Ref. 12(| . The latter feature is, according to Voitkiv et al., a clear signature of the 
ee- Auger mechanism. Thus, our numerical calculations that account for all the first-order 
(nonradiative) mechanisms using the accurate, highly correlated wave function of helium do 
not support the findings of Refs. 10Ml2|. In this connection, it should be noted that angular 
correlations due to the interaction V12 play very important role, if being included in the trial 
helium wave function $(r*i,r2). They strongly influence the momentum distribution of the 
emitted electron in the backward direction, which is due to the shake-off mechanism 21 ]. 
And in Fig. 1 we can see manifestation of their effect in the region k z < 0. However, it is 
quite different from manifestation of the ee- Auger mechanism claimed by Voitkiv et al. It is 
interesting that the ee- Auger mechanism has no manifestation in the OBK approximation 
either, as can be concluded from Fig. 2. According to Eq. (9), the effect of Vyi found in 
Refs. 10Nl2[| is clearly cancelled by the other first-order mechanism, which involves the target 
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nucleus (the Vjvi term). This finding markedly illustrates the importance of accounting for 
all the binary interactions between the particles taking part in the reaction. 

V. SUMMARY AND CONCLUSIONS 

We considered theoretically a transfer ionization channel in a fast proton-helium collision, 
focusing on the so-called first-order electron-electron capture mechanism proposed recently 



by Voitkiv et al. |10Hl2|. It was shown, using consistent quantum scattering theory, that this 
specific mechanism is included in a well known OBK scenario. The first-order calculations 
with an accurate trial helium function exhibited no signature of any significant effect of the 
ee- Auger process. 

Some comments should be made with regard to the above theoretical analysis. The 
equivalence of the post- and prior-forms of the transition amplitude is realized only if the 
exact helium wave function is employed. In that case, we have the following equation for 
this function: 

(e% e - h w - h 20 - V N2 )\%) = (Vm + V 12 )\%). (17) 

It can be readily shown, using energy conservation and properties of the final asymptotic 
state, that the projectile-electron potential V p i in the left-hand side of Eq. (9) can be replaced 
with the operator (e^ e — hio — ^20 — Vm) occurring in the left-hand side of Eq. (17). Thus, 
the left- and right-hand sides of the Schrodinger equation (17) correspond to, respectively, 
the prior and post matrix elements in Eq. (9). This feature explains a well documented 
fact that Eq. (9) is fulfilled to a good approximation in the case of an accurate trial helium 
function, which is typically obtained from a variational procedure. But if the trial function 
is poor, then the right-hand side of Eq. (17), which is related to the post matrix element 
in Eq. (9), appears to yield a better approximation to the exact result than in the case of 
the left-hand side. This observation explains why the results of Voitkiv et al. using the 
correlated and uncorrelated helium functions in their post-amplitude are similar (see Fig. 7 
of Ref. m 



Finally, in their calculations Voitkiv et al. 1CH12J] adopt an approach in which they 
incoherently sum contributions to DDCS from different mechanisms, including, in particular, 
the OBK (or capture- shakeoff) and ee-Auger. However, such an approach is manifestly 

11 



incorrect because, as follows from the above theoretical analysis, the ee-Auger contribution 
is already included in the OBK approximation. Thus, the corresponding numerical results 



of Refs. 



10h12| can be regarded as erroneous. 
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Figure 1: FBA distribution of the electron momentum components at E p = 3.6 MeV. DDCS is 
given in barn/(a.u.) 2 
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Figure 2: The same as in Fig. 1, but in the OBK approximation. 
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